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Abstract. Internet is known to display a highly heterogeneous structure and complex fluctuations in its 
traffic dynamics. Congestion seems to be an inevitable result of user's behavior coupled to the network 
dynamics and it effects should be minimized by choosing appropriate routing strategies. But what are the 
requirements of routing depth in order to optimize the traffic flow? In this paper we analyse the behavior of 
Internet traffic with a topologically realistic spatial structure as described in a previous study (S-H. Yook 
et al. ,Proc. Natl Acad. Sci. USA, 99 (2002) 13382). The model involves self-regulation of packet generation 
and different levels of routing depth. It is shown that it reproduces the relevant key, statistical features 
of Internet's traffic. Moreover, we also report the existence of a critical path horizon defining a transition 
from low-efficient traffic to highly efficient flow. This transition is actually a direct consequence of the 
web's small world architecture exploited by the routing algorithm. Once routing tables reach the network 
diameter, the traffic experiences a sudden transition from a low-efficient to a highly-efficient behavior. It 
is conjectured that routing policies might have spontaneously reached such a compromise in a distributed 
manner. Internet would thus be operating close to such critical path horizon. 

PACS. 89.75.-k Complex Systems - 05.70. Ln Nonequilibrium and irreversible thermodynamics - 87.23. Ge 
Dynamics of social Systems 



1 Introduction 

The efficient performance of any communication network 
, is jeopardized by congestion problems, which often show 

■ up in unpredictable ways. This seems the case, for exam- 
" pie, of so called Internet storms Such problems were 

early identified in different types of engineered networks. 
Norbert Wiener for example mentions that "a switching 
, service involving many stages and designed for a certain 

■ level of failure shows no obvious signs of failure until the 
' traffic comes up to the edge of the critical point, when it 

goes completely into pieces, and we have a catastrophic 
traffic jam" j^. These observations allow to formulate a 
, number of key questions concerning communication nets, 
' such as: How do critical traffic levels are reached? What 
is the nature of these thresholds? How appropriate rout- 
ing algorithms modify this behavior? Are there optimal 
routing strategies? 

Modeling Internet dynamics has been an active area 
over the last decade. The approaches include detailed sim- 
ulations [3, simple statistical models and verbal models 
In this context, in |S] we investigated the existence of a 
jamming phase transition between the free phase and the 
congested phase in a model of network traffic over regular 
meshes. This phase transition depends on traffic density. 
It was conjectured that large-scale network traffic self- 
organizes near the critical point of the transition, which 
is linked to high network efficiency and unpredictability 




Fig. 1. Small Internet-like network (see text for generation 
method description). The parameters have been set to the po- 
sition corresponding to Internet in the (D/, a, a)- phase space: 
M= 250, < fc >= 4, _D/ = 1.5, Q = 1, cr = 1. 



At the critical regime, the distribution of congestion 
duration lengths scales as a power-law j2j jZi ■ 

Such a self-organized scenario is reinforced by a recent 
study suggesting that Internet fluctuations are a conse- 
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quence of the internal dynamics of the system . Against 
previous claims, there is mounting empirical support that 
network traffic heterogeneity is a consequence of collective 
dynamics and not because of the high variability injected 
by external sources. 

In this paper we will show that Internet is able to route 
efficiently its inner flow of packets because a special com- 
bination of local routing rules and a particular network 
architecture. As we will see, Internet routing reaches an 
optimal, low-cost traffic flow as a result of a trade-ofF be- 
tween random and deterministic routing schemes. 



2 Modelling Internet's traffic 

Following previous approaches |n]|Z| let us consider net- 
work defined on a graph fl with M nodes (figure 1). The 
number of links connecting a given node with another (ir- 
respective of their characteristics) will be indicated as ki . 
The subset of closest nodes oi Si ^ Q is Q. Only a fraction 
of fixed p < (0, 1] nodes is selected as source and desti- 
nation traffic endpoints (hosts) . Host locations are chosen 
randomly. The other (1 — p)M nodes can only store and 
forward incoming messages (routers). 

Both host and routers can route only one packet at a 
time (the routing policy is described below). Each node Si 
is provided with a finite queue of packets waiting for com- 
munication resources (free links). Queues are not allowed 
to store more than H packets simultaneously. New pack- 
ets arriving at an already saturated queue will be simply 
removed from the system. This provides a very simple dis- 
sipative mechanism useful for recovering from heavy traffic 
congestion. If n{si,t) is the number of packets at Si, the 
total number of packets in the system will be 

7V(i)= ^ n(s„i) (1) 

which is the key quantity which has to be analysed here. 

Any microscopic host behavior compatible with the 
fluid scenario must ensure a controlled injection of pack- 
ets in order to not flooding the network (and thus entering 
the congested state). Note that in a fluid traffic regime it 
is unlikely that packets will be alive for an arbitrarily long 
time. A fluidity requirement will necessarily impose a hard 
constraint on the maximum time spent by a packet trav- 
elling across the network. The simplest way of ensuring 
fluid packet flow is to stop the sources emitting new pack- 
ets when detecting local congestion, that is, when there is 
no empty space in the source neighborhood C'i to fill with 
new packets. If no self-regulation is present, it has been 
shown that there is a sudden and sharp jamming transi- 
tion from the free to the congested state depending on the 
mean (system) packet generation rate (A) . At the critical 
point between the two phases the system displays opti- 
mum (global) performance. If self-regulation is allowed, it 
can be shown |7| that simple traffic source rules are able 
to self-organize the traffic network around a definite mean 
rate (A), which scales as the quotient between the average 
system load < N > and the averaged packet latency (T) : 



where the latency is defined as the time comprised from 
the creation of a packet until its delivering at the destina- 
tion host. The mean latency < T > is averaged over all 
successfully released packets. Following [7] we will indicate 
by ^ the number of local congested neighbors: 

e= ^K^i)] (3) 

/c6C(i) 

where 9[x] = 1 for a; > and zero otherwise. The packet 
injection rate for host i is updated as follows: 

{ mm {1, Xi{t) + n} ^ = 
Kit) Q<^<h (4) 

^^k, 

where fi is the so-called driving parameter. The sec- 
ond rule allows a particular host to stabilize around a 
given rate. Traffic rate increases conservatively and drops 
down to zero when all neighboring nodes are congested. 
The reader must be aware that the above rules are not in- 
tended to be a detailed model of real traffic sources. Traffic 
sources can not be described with an universal distribu- 
tion probability. Moreover, it seems that a rich variety of 
distributions (exponential, bi-modal or log-normal) apply. 
As we will see later in the paper, it is unlikely that the 
variability of traffic sources will be the cause of the scal- 
ing detected in Internet traffic. Moreover, this model does 
not introduce any explicit correlations between sources 
(i.e.: like in active conversation between two hosts) so 
the dynamics can not be the by-product of pre-defined 
source correlations. Because the system is poised to crit- 
icality (the so-called fluid packet flow regime) the global 
dynamics emerge from the collective behavior of its com- 
ponents. Ultimate source details are irrelevant in this con- 
text. Moreover, although TCP is more realistic than our 
simplifled local rules (also much more complex) our con- 
jecture is that the large-scale dynamics will not be greatly 
changed by the local behavior of hosts within this fluid 
flow regime. This idea has been partly answered in a re- 
cent work in which it has been shown that IP level 
traffic (packet dynamics) does not depend to any signif- 
icant extent on the TCP arrival process (host behavior), 
supporting our view that traffic dynamics is not a conse- 
quence of detailed source (host) behavior. 

For a homogeneous network with average degree < k > 
(and finite < fc^ >) it can be shown, following [71 that the 
time evolution of packet density r{t) — N{t)/M in the 
limit of _ff ^ cxD is defined by the mean field equation: 

-.,._lA^ra-r) (5, 

where D indicates the network diameter (also the average 
transient time, assuming packets jump from node to node 
with the same average time scale). It is not difficult to 
show that the equilibrium points of the previous equation 
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are given by r± = [1 ± (1 - Ap\D / < k >y^^]/2. For 
A > Ac =< k > /4:Dp, the fixed points vanish and no 
finite density exists. In this "congested phase" the density 
of packets grows without bounds. For A < Ac, a finite 
stable density 
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is observable (the other fixed point r+ is unstable). For 
the particular case p — 1 analysed by Fuks and Lawnizak 
jl2|. we recover their critical point Ac = 2/L for dynam- 
ics taking place on a square lattice, i. e. D — L/2. If a 
feedback exists between A and 7, some finite equilibrium 
density F* will be achieved in the previous model and thus 
no divergence will be allowed to occur. In this case, the 
mean field model indicates that a scaling relation will be 
observed, i. e. 

A^p-i (7) 

between the (self-organized) packet release rate and the 
density of hosts (which is here the only relevant external 
parameter). Such a scaling relation was shown to occur in 
the previous model. 

The previous mean field calculation was done by as- 
suming that a fixed A is beign used. The previous rules ac- 
tually introduce self-regulation of injection rates by traffic. 
In other words, if traffic level N defines an order parame- 
ter, it will interact with a control parameter (A), reducing 
it when N is large and incresing it when low. The feed- 
back between these two key quantities results in a state 
dominated by fluid, but fluctuating traffic with many char- 
acteristics in common with observed Internet's dynamical 
patterns. 

In order to explictly consider the feedback between or- 
der and control parameters, we can consider a new mean 
fleld approximation based on the previous local rules. It 
can be shown, asuming finite H , that the new set of equa- 
tions is now: 
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for low density levels (i. e. _r <C i?, consistently with a 
fluid traffic) the single fixed point (obtained from dF/dt — 
dX/dt = 0) is 
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the Jacobi matrix L for the previous set of equations is 
given by 

P 



-<k> 
D 
-1 
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(11) 



The associated eigenvalues are 
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Fig. 2. Network dynamics and routing: it involves a given 
depth of routing m (a path horizon). A packet traveling from Sj 
to Si within the 5'"(si) domain (of depth m) is deterministically 
routed along the shortest path {d{j,i) < m). The packet tra- 
verses hosts (squares) and routers (circles) indistinctly. Packets 
traveling outside the m-domain (i. e. for d{j, i) > m) have more 
than one path choice and perform random walks. As soon as 
the packet enters into the m-domain, the packet is determin- 
istically routed along the shortest path. Here we would have 
m = 2. 



both of them real and negative: the point attractor is glob- 
ally stable. Numerical simulations of the model on a Pois- 
sonian graph (where the previous approximation would 
held) agree with these predicted values. Topological fea- 
tures are included only as averaged quantities, here mean 
degree < fc > and diameter D. However, our interest is 
to explore the traffic dynamics on a realistic network ar- 
chitecture, in order to provide the closest modeling ap- 
proach under the previous rules. It has been shown that 
the Internet does not display the homogeneous architec- 
ture assumed by the Poissonian graphs '13' . In the next 
section the network's topology used in our analysis is pre- 
sented. Together with an explicit definition of the routing 
algorithm, they will complete our model's description. 



3 Network topology 

In previous analyses [H] [Zl the topology of the graph was 
chosen to be a lattice. Although this might seem a limited 
topological arrangement, it has been successfully tested 
on real hardware and the presence of a phase transition 
fully confirmed 14 . Moreover, mesh connected topologies 
may be the most efficient solution at the limit of very large 
parallel computer architectures dj ■ 

Regular lattices fail to describe Internet topology. Re- 
markably, Internet displays a scale- free architecture 
The origins and causes of this topology have been the sub- 
ject of most discussion and controversy. It has been show 
that most existing Internet generators fall in a very differ- 
ent region of the phase space where real Internet is located 

This phase space is defined by (-D/, cr, a), where every 
parameter defines a major force shaping a different as- 
pect of the large-scale Internet topology. Note that this 
model does not define all detailed correlations observed in 
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Internet and/or the precise functional form of Internet's 
path length and degree distribution (i.e.: the exact expo- 
nent of the scale- free distribution) . This is a minimal set 
of universal parameters that any realistic Internet model 
must satisfy and it is a very good approximation to the 
large-scale topology. 

The spatial distribution of Internet nodes is not ran- 
dom. It has been noticed a strong correlation between 
fractal distribution of cities and Internet nodes. The mea- 
sured fractal dimension is Df = 1.5. When generating the 
Internet-like network, the position of nodes is obtained by 
sampling a Rayleigh-Levy dust of the same fractal dimen- 
sion. The likelihood of placing a link between two nodes 
Si and Sj depends both on the (euclidean) link length Wij 
and linear preferential attachment: 

7T(fc,,d,,)«-^ 

Note that longer links cost more and thus will be se- 
lected with less probability. Traditional topology gener- 
ators based on Waxman model J7| wrongly assume expo- 
nential decay instead of linear cost decay. On the other 
side, there is empirical evidence that highly connected 
nodes will be linked with higher probability. Increasing 
a will favor linking to nodes with higher degree, while a 
higher a will penalize longer links. From real measure- 
ments PP , it has been identified the position of Internet at 
Df = 1.5, a = 1 and a = 1. Here, all numerical simula- 
tions have been performed using this topological arrange- 
ment (figure 1). 

4 Path horizon and routing tables 

An important ingredient when modeling network dynam- 
ics is to describe the paths followed by packets towards 
their destinations, that is, the routing policy. By prop- 
erly defining this routing algorithm, we will complete our 
model's definition. 

Real routing protocols do not drive packets at random. 
Instead, they try to route packets along the most efficient 
routes (i.e: minimize distance or latency). At the same 
time, it is unrealistic to assume that all packets follow 
optimal paths because the large amount of global infor- 
mation replicated at every single node. Clearly, the paths 
traced by packets in real networks are properly character- 
ized as a trade-off between random diffusion and optimal 
routing. This is reflected in the real Internet by its two- 
level routing organization. Nodes are grouped in so-called 
Autonomous Systems (AS) and different routing rules ap- 
ply at each level. Intra-AS routing is based on shortest 
path routing but inter-AS routing does not cleary follow 
any minimization criteria. 

A simple way to explore the cost/efflciency trade-off is 
by introducing a parameter defining the visibility scope of 
the node (depth of routing parameter m or node domain 
diameter), that is, a sphere /^(™^(si) of radius m centered 
at every node Si G H (figure 1). We allow every node to 
known every other node at a distance of m hops or less 




Fig. 3. Hierarchy of path sets defined by the routing pohcy. 
Every routing scheme explores a fraction of network paths (a 
subset of the entire path set P(Q)). Pure random walk strat- 
egy (m=0) visits aU available paths, even shortest paths. In 
fact, random walk travel more frequently along those shortests 
paths. The more restrictive set is associated to pure shortest 
path routing (m > [(d)]), which chooses a small fraction of all 
possible routes on the network. 



but no more. No information will be stored about nodes 
outside the node domain. This idea is indicated in figure 
2, where a given target node Si is shown at the center 
of its TO-sphere. When a foreign node Sj & H — /^(™'(si) 
send a packet towards Si, while moving in the outside of 
the sphere it performs as a random walk. Once the packet 
hits the boundary of the sphere, dr'-"^'> (si), it is routed 
along the shortest path JH] • 

5 Efficiency and network's exploration 

The depth of routing parameter m induces a hierarchy 
of path subsets over the entire set of available network 
paths (see figure 3). The random m = and deterministic 
m — M routing represent the most general and restric- 
tive subsets, respectively. Increasing m will progressively 
reduce the randomness at routing decisions. Here, we find 
a nested collection of subsets corresponding to the inter- 
mediate situations < m < M. The relative size of each 
subset can be approximated by measuring the fraction of 
nodes F visited by packets: 

1 " 

^=mE^[^'] (12) 

i=l 

where Ai is the probability of visiting the node Si. This 
term depends both on the network dynamics and the ge- 
ometrical situation of the node within the network. Con- 
siderable effort has been devoted to characterize the like- 
lihood of visiting a node in a purely static fashion. In this 
context, two relevant centrality (or 'load' |20j(2T]) mea- 
sures are Random Walk Centrality [22] and Betweenness 
Centrality '23' which certainly represent the two extremes 
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4 6 

Path horizon m 

Fig. 4. The fraction of visited nodes depends on the amount 
of order defined by the depth of routing parameter. Determin- 
ism sharply constraints the routing paths. Every point was 
obtained averaging over an ensemble of four different networks 
and every single network was simulated in five different host 
configurations. Network parameters: M = 512, < k >= 4, 
p = 0.1, II = 0.01, H = 10, Df = 1.5, cr = 1, a = 1 and 
r = 7 X 10^ steps. 



of our hierarchy of path subsets. Anyway, the correlation 
of these measurements with dynamic centrality Ai is weak, 
to say the least. 

In figure 4 the numerically measured F{m) for a fi- 
nite range of m values is shown. This fraction equals 1 
for random routing < to < [{d)], where (d) is the av- 
eraged shortest path length and [x] denotes the integer 
part of X. Random routing does not discard any network 
node. A sharp transition takes place from this point [{d)] < 
m < M and a large fraction of network nodes (about 45 
percent) are never visited by deterministic routing. This 
severly restricts the diversity of routes and yields a load- 
insensitive system |24) . The system requires a certain de- 
gree of noise in order to avoid sending packets through 
already collapsed nodes while it is enabled to choose less 
optimal but free routes. Note that depth of routing pa- 
rameter TO defines an order parameter because quantifies 
the degree of order existing in the system. The existence of 
an order parameter confirms the critical nature of fnternet 
traffic. 

Numerical simulations have shown that flow is max- 
imized at the order-disorder transition point m — [{d)]. 
This can be observed in network throughput reaching a 
maximum at this point (see figure 5). Throughput is de- 
fined by the quotient of the number of successfully released 
packets and the sum of all packets generated during the 
simulation. In the next section, we will show that several 
real statistics can be reproduced by the system dynamics 
poised to criticality. 




Fig. 5. Exploring the network traffic dependency on depth 
of routing parameter (m). The lines connect the points for il- 
lustrating purposes only. (A) Mean latency is considerably re- 
duced when enough system information is given to individual 
nodes. Most packets are forwarded along the minimum num- 
ber of hops. (B) Global throughput is optimal at the criti- 
cal point when path horizon is about the network diameter. 

(C) Mean packet rate also drops down at intermediate value. 

(D) Mean workload also experiences sudden transition from 
heavy to light load. Network parameters: M=250, < k >= 4, 
Df — 1.5, (T = 1, Of = 1. Simulation parameters: p = 0.1, 
fi = 0.01, H = 10, T = 10^ steps. The shape of these plots 
does not depend on network size, that is, the optimal point is 
always found at m = [{d)]. 



6 Average stretch 

The average stretch s measures the efficiency of routing 
by comparing the number of hops h traversed by a packet 
to the shortest path distance d between source and desti- 
nation: 



s = -; 



Compact routing schemes minimize the average stretch 
while maintaining the size of routing tables small I25J • Re- 
ducing the average strech will progressively raise up the 
memory requirements at each router. Assuming Thorup- 
Zwick (TZ) compact routing scheme jlHI: the average stretch 
can be expressed as a function of the distance distribution 
and the graph size M only: (s) — f{{d) ,ad) |2IZ|- This 
TZ scheme ensures a nearly optimal lower memory upper 
bound for (s) = 3 in generic networks. For scale-free net- 
works, TZ achieves lower bounds. In particular, for the 
Internet interdomain graph (Autonomous Systems) with 
degree distribution exponent ^ k, 2.1 and M = 10'*, the 
TZ average stretch (s) « 1.14p7j. Also, the average num- 
ber of entries in the routing tables is approximately 52. It 
turns out that (s) {{d) , Ud) surface has unique minimums. 
Strikingly, the points corresponding to Internet distance 
distribution are very close to them I27| . This suggests that 
Internet topology is shaped by some hidden optimization 
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Fig. 6. (A) Average stretch (B) Average fraction of nodes 
inside m-domain is an estimation of the amount of informa- 
tion required by each node. Here open circles correspond to 
the Internet's nodel. For comparison, the simulation has been 
repeated in a poissonian network of same size and mean de- 
gree (filled circles). Network parameters: M = 500, < k >= 4, 
Df — 1.5, a = 1, a — 1. Simulation parameters: p — 0.1, 
/i = 0.01, H = 10, T = 5 X 10* timesteps. Measurement window 
= 200. The distributions were obtained from an statistical en- 
semble of three networks, every network simulated three times 
with hosts located at different configurations. 



criteria. Anyway, TZ scheme is not a realistic Internet in- 
terdomain routing scheme because assumes that global 
topology view is available. We have numerically measured 
the average stretch and the average fraction of neighbours 
at distance m for our routing scheme (see figure 6). Note 
that the critical path horizon m — [{d)] is very close to 
the minimum average stretch. 



7 Network fluctuations and performance 

In order to test the goodness of the model presented here, 
it is interesting to compare some real Internet statistics 
with their respective model measurements. In particular, 
those measurements will be collected for the model at the 
critical path horizon, that is, when m = [{d)]. 

When understanding the competition between a net- 
work's internal collective dynamics (i.e: Internet traffic) 
and external environmental changes (i.e: traffic sources or 
host behaviour), it is useful to study the relationship be- 
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Fig. 7. The relationship between fluctuations and the average 
incoming node flux (a similar distribution holds for average 
outgoing router flux). The plot shows that both quantities are 
related by a power law of exponent 1/2, which is consistent 
with the measurements from Internet routers. Network param- 
eters: M=500, < k >= 4, Df = 1.5, a = 1, q = 1. Simulation 
parameters: p = 0.1, p. = 0.01, H = 10, T = 10000 steps. Mea- 
surement window = 200. The distributions were obtained from 
an statistical ensemble of 10 networks, every network simulated 
5 times with hosts located at 5 different configurations. 



tween the mean flux and the size of fluctuations around the 
average Previous explanations of ocurrence of self- 
similarity in traffic networks are based on the superposi- 
tion of many and high- variable (infinite variance) sources 
This point of view discards the effects of the sys- 
tem's collective dynamics and considers that Internet is 
an externally driven system. Real data shows that Inter- 
net dynamics can not be simply reduced to the behaviour 
of traffic sources. 

Let us define the incoming (outgoing) flux fi as the 
amount of packets received (forwarded) at router Si during 
a given and fixed period of time. For every router, compare 
the average fiux (fi) with the dispersion cr^ around the 
mean. It has been noticed that for several real systems 
the following scaling relation holds: 



if)' 



where a is an exponent which can take the values of 
1/2 and 1 jTHj. This suggests that real systems can be 
classified in two main classes depending on the value of 
this exponent. The relevant exponent in this context is 
a = 1/2, which has been observed in daily traffic mea- 
surements at 374 geographically distinct Internet routers 
|l()j . Systems exhibiting the 1 /2 exponent are representa- 
tive of endogeneous dynamics, that is, determined by the 
system's internal collective fluctuations. Moreover, mea- 
surements on our model reproduce the same exponent (see 
flgure 7). 

A measure of Internet end-to-end performance is the 
normalized latency time r |28j deflned as the quotient be- 
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tween latency time L (measured as Round- Trip-Time) and 
geographical distance w: 

_ L 
w 

There are several factors governing r. First, propa- 
gation rate is finite. There is a minimum delay because 
packets can not move faster than speed of light. Second, 
the number of nodes traversed departs from the mini- 
mum found along the geodesic path from node to destina- 
tion. And third, packets spent some time at every inter- 
mmediate node because queueing delays. Define Tmin — 
Lmin/ui as the normalized latency without taking into 
acount queue delays and Tq„ — Lav/w as the normalized 
latency while considering all factors affecting packet la- 
tency. 

Their probability distributions have been measured from 
two years of PingER [221 data and follow power-law scal- 
ing with stable exponents of about -3.0 for and -2.5 
for Tav p8_. Numerical simulations on the model poised 
to criticality reproduced the power-law P{Tav) probabil- 
ity distribution. The exponent of the distribution is about 
-2.45, which is very close to real observations (see figure 
8). The previous results seems to be quite robust and in- 
dependent of most model parameters changes. An inter- 
esting thing to note is that current model cannot repro- 
duce the existing correlation between Tmin and Tav and 
geographical distance w reported by previous studies f2Bi • 
The reason might be that propagation rate is (unrealis- 
tically) assumed to be infinite in the model. It might be 
that this factor has little influence in shaping the Tav prob- 
ability distribution. Moreover, it is worth noting that any 
deviation of the order parameter from the critical regime 
resulted in an exponential distribution for Tav, reinforc- 
ing the view of an optimally efficient and self-organized 
Internet. 



8 Discussion 

In this paper we have shown that a simple model of traf- 
fic dynamics incorporating the appropriate Internet's net- 
work topology is able to recover several statistical fea- 
tures of real traffic. More important, we have seen that the 
routing algorithms can take advantage of the small- world 
structure of the web by reaching a critical path horizon 
close to the network's average path length. In doing so, 
a highly efficient system is reached at low cost: routing 
strategies only need to consider a small depth. Once the 
m parameter reaches the network's diameter, no further 
information is required to properly reach the target. A 
full-system deterministic routing strategy is actually un- 
necesary and would be too costly. Instead, the constraints 
imposed by network's architecture allow to exploit the im- 
plicit information defined by the small-world architecture. 

It might be that Internet evolves in a way that through- 
put or global performance is maximized. This trend is con- 
strained within the limits of available communication re- 
sources. Inside this regime, Internet is shaped in order to 
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Fig. 8. Normalized latency Tav distributions at the critical 
point m = [(d)] follows power law of exponent ~ —2.45 and 
fits very well the real measurements (see text). Here, we plot 
the distributions for H — 5 (shaded line) and H — 10 (contin- 
uous line) showing that the long tail does not depend on the 
maximum queue size. Any deviation from m = (d) results in 
an exponential distribution, deviating from real observations. 
Network parameters: M = 1000, < k >= A, Df = 1.5, a = 1, 
a = 1. Simulation parameters: p — 0.1, /i = 0.01, H = 5,10 
and r = 5 X 10^ steps. The distributions were obtained from an 
statistical ensemble of 10 difi'erent networks and every network 
was simulated five times with different host arragements. 



provide better response. This is reflected in the optimal 
and lower average stretch observed in real Internet, which 
is close to global minimum j^Hl- This suggest hat some 
hidden optimization is at work. Clearly, an unresponsive 
system will be no good. How a distributed collection of 
designers were able to define this globally efficient infras- 
tructure is a question that deserves attention. 

This paper is dedicated to the memory of Per Bak. This work 
was supported by a grant BFM2001-2154 and by the Santa Fe 
Institute. 
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